arXiv:1506.00049v2 [astro-ph.HE] 16 Dec 2015 


Two-stage acceleration of interstellar ions driven by high-energy lepton plasma flows 

cm YunQiani, SHENG ZhengMing^-S-^j] LU QuanMing^, LI YuTong^-^, and ZHANG Jie^d 

^ Beijing National Laboratory of Condensed Matter Physies, 

Institute of Physics, CAS, Beijing 100190, China; 

^ SUPA, Department of Physics, University of Strathclyde, Clasgow Cf ONC, United Kingdom; 

^ Key Laboratory for Laser Plasmas (MoE) and Department of Physics and Astronomy, 

Shanghai Jiao Tong University, Shanghai 200240, China; 

^ IPSA Collaborative Innovation Center, Shanghai Jiao Tong University, Shanghai 200240, China and 
^ School of Earth and Space Sciences, University of Science and Technology of China, Hefei, 230026, China 

(Dated: December 17, 2015) 

We present the particle-in-cell (PIC) simulation results of the interaction of a high-energy lepton 
plasma flow with background electron-proton plasma and focus on the acceleration processes of the 
protons. It is found that the acceleration follows a two-stage process. In the first stage, protons are 
significantly accelerated transversely (perpendicular to the lepton flow) by the turbulent magnetic 
field ’’islands” generated via the strong Weibel-type instabilities. The accelerated protons shows a 
perfect inverse-power energy spectrum. As the interaction continues, a shockwave structure forms 
and the protons in front of the shockwave are reflected at twice of the shock speed, resulting in a 
quasi-monoenergetic peak located near 200MeV under the simulation parameters. The presented 
scenario of ion acceleration may be relevant to cosmic-ray generation in some astrophysical environ¬ 
ments. 
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INTRODUCTION 


The origin of high-energy cosmic-rays over a wide 
range of energies and the non-thermal emission of ra¬ 
diation from a wide variety of high energy astrophysical 
sources is the most fundamental problem in astrophysics 
and has been studied for over six decades since Teller 
and Fermi [l[, however, it has not been fully resolved till 
now. Usually these radiations are attributed to mecha¬ 
nisms like inverse Compton scattering 0, Q and/or syn¬ 
chrotron emission 0 from high energy particles and the 
problem turns to how these power-spectrum energetic 
particles are generated. Considerable efforts are devoted 
to finding possible generation mechanisms of such ener¬ 
getic particles [1, As far back as 1949, Fermi pur¬ 

posed the stochastic acceleration mechanism suggesting 
that the particles can be accelerated through the colli¬ 
sion with magnetic ’’islands” in the space [l|. Recently, 
Hoshino presented a similar process using a magnetic re¬ 
connection configuration in pair plasmas Q- Another 
hopeful candidate of ” cosmic particle accelerator” is the 
collisionless shock. It is believed to exist widely in the 
interstellar space and play an important role in super¬ 
nova remnants (SNRsiM l^. jets of radio galaxies 
gamma-ray bursts (GRBs)[ll|, ll2| and the formation of 
the large scale structure of the Universe [Tsi Q- There¬ 
fore, it has been modelled by pl enty of researchers both 
analytically [l5| and numerically [ igI fl^ . 

Although both stochastic acceleration and collisionless 
shock acceleration have been well modelled separately 
in lots of publications, most of them are restricted to 


the phenomena themselves. Presumed particular initial 
conditions are used to ensure the occurrence of the in¬ 
terested phenomenon. For example, in Hoshino’s work 
0 , four Harris current sheets are imposed to build niag- 
netic islands; in the collisionless shock simulations [16- 
El, a reflecting boundary condition is used to ’’make” 
the shock. These studies have significantly improved our 
general understanding of the stochastic acceleration and 
shock dynamics. However, more realistic modelling in a 
wide parameter range is still needed. 

In this work, we present a scenario of proton ac¬ 
celeration processes in two distinct processes occurring 
naturally during the interaction of a high-energy lep¬ 
ton plasma flow or jet with background electron-proton 
plasma. This system of interaction corresponds to sev¬ 
eral astrophysical models such as the pulsar-wind shock 
systemfioj. the fireball model of GRBsjiH and the re¬ 
cently pi^osed Binary orbits model of novae y-ray 
emission |22|. Also the electron-positron jets are found 
in different astrophysical environments such as quasars 
[ 2 ^ , black holes , and so on. Without using any arti¬ 
ficial setup to facilitate any specific acceleration process, 
by use of the particle-in-cell (PIC) simulation, we find 
that the background protons in the system are naturally 
accelerated via Fermi acceleration like processes. The ac¬ 
celeration can be divided into two stages which are dom¬ 
inated by different mechanisms [25[. In the early stage, 
stochastic acceleration, which is driven by turbulent mag¬ 
netic fields generated by strong Weibel-type instabilities, 
dominates and the protons are accelerated transversely. 
As time elapses, a collisionless shock gradually forms in 
the background plasma, which accelerates the protons 
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longitudinally. These results provide a clear physics pic¬ 
ture of a high-energy plasma flow (jet, pulsar-wind, and 
so on) interaction with interstellar masses and thus can 
improve our understanding of the related astrophysical 
phenomena. 


SIMULATION SETUP 

The numerical simulation was performed by our self- 
encoded two-dimensional PIC code KLAPS [1^. We 
simulate a two-dimensional unmagnetized system with 
the electron-ion background plasma which homoge¬ 
neously fills the whole simulation domain located at 
{x,y) G [0,900(ie] x [0,95(ie], where dg = c/ujpe = 
{ATTe^nplmeC^)~^^‘^ is the election skin depth for a num¬ 
ber density Up and electron mass rrie] e and c here de¬ 
note the elementary charge and the speed of light in 
vacuum, respectively. We use a real ion-electron mass 
ratio in the simulation, i.e., mijm^ = 1836. The ini¬ 
tial temperature for background plasma is 500eV. A mo- 
noenergetic pair plasma flow consists of electrons and 
positrons which are injected along direction from the 
left boundary of the simulation domain with initial en¬ 
ergy ^ 50MeV(7e+g- = 100) and interact with the back¬ 
ground plasma. The initial density of the flow n/ is equal 
to Tip. The whole system is numerically resolved with 8 
cells per dg in both directions with 9 particles per cell for 
each species, thus ensuring that the behavior of leptons 
can be modelled accurately. The time resolution of the 
simulation is 0.06uj~J-. An open boundary condition is 
used for x-direction and a periodic boundary condition 
for ^-direction. 


RESULTS AND DISCUSSIONS 
First stage: stochastic acceleration 




FIG. 1. The space distributions of (a)electrons, (b)positrons, 
(c)ions and (d)magnetic fields as well as {e)x — px, {f)x — py 
and {g)px —Py projections of the ion’s phasespace distribution 
at t = 950ct;“e^. (h) shows the energy spectrum of background 
ions (blue line) and electrons (orange line) at the same time, 
with the line ~ (dashed line). 


celerated in both x and y directions. The longitudinal 
acceleration in x direction can be simply attributed to 
the momentum transfer between the flow leptons and the 
background particles through the electromagnetic fields. 
However, it can be seen from Figure [TJg) that the trans¬ 
verse acceleration is even stronger than the longitudinal 
acceleration. Moreover, the accelerated protons forms a 
perfect inverse-power energy spectrum 

dN/dE - E-^ (1) 


In the beginning, the flow is uniform along the y- 
direction. However, due to the relative motion between 
the lepton flow and background, Weibel-type instabil¬ 
ities can occur and generate filamentous structures at 
X < lOOdg after several hundreds of l/cjpg. These fila¬ 
ments make the newly incoming flow leptons filamented, 
as shown in Figure [T](a)-(b). This filamented high-energy 
lepton plasma flow then propagates deeper and causes a 
highly turbulent region at 100-400dg after about 900co’“g^. 
Within the turbulent region, the density of ions can vary 
more than three orders of the magnitude, as shown in 
Figure [TJc). Such strong turbulence can generate strong 
magnetic fields which are displayed in Panel (d) of Figure 
[1] [27|. As said before, the paper mainly focuses on the 
acceleration of the ions in the background plasma (re¬ 
ferred as ’’ions” after for simple). Figures [TJe) -(g) gives 
the phasespace distribution of the ions. The ions are ac- 


with the cutoff energy near 100 MeV, as shown in Figure 

HTh). 

To interpret the strong transverse acceleration and 
the inverse-power energy spectrum, we notice the strong 
turbulence-generated magnetic fields. Fermi demon¬ 
strates that particles can be accelerated by collisions with 
the moving magnetic islands, and proves that for a large 
number of collisions, the resulted energy spectrum will be 
in the form of dN/dE ^ E] where r is the 

average time between two collisions and T is the lifetime 
of the particle before being absorbed or escaped from the 
region. This is often referred to as Fermi H acceleration. 
Consider the following fact that (1) our code for this sim¬ 
ulation does not include the annihilate or recombination 
module, i.e., the local particle number conserves, (2) a 
periodic boundary condition is used in ^-direction, and 
(3) the x-direction size of the simulation domain is much 
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FIG. 2. The momentum space trajectory of 100 ions from 
background plasma with the largest transverse momentum at 
t = 950a; 


larger for ions initially located at x < dOOde to escape. 
It is thus reasonable to assume that T ^ oc and then 
dNjdE ^ which is exactly the result we get in the 
simulation. We also trace the trajectories of the 100 ions 
in momentum space with the largest transverse momen¬ 
tum at t = 950a;as plotted in Figure [2l One can 
see obvious evidence of the reflection between significant 
acceleration or deceleration marked by the grey circles, 
which indicate that some ions experience sequential ac¬ 
celeration occurs at different locations. This is different 
from the usual acceleration mechanisms of ions found in 
relativistic laser-solid interaction in laboratory [28| . 

We note that some researchers have reported some 
results from similar models with much smaller incident 
flow energy {'y/iow < 10) since one decade ago and have 



FIG. 3. (a) The space distribution of ion density, (b) the 

transverse averaged ion density, the space distribution of (c) 
magnetic field Bz and (d) electric field Ex, (e) transverse aver¬ 
aged field energy densities, blue for magnetic field Bz and red 
for electric field Ex, {f)x—px and {g)px —Py projections of the 
ions’ phasespace distribution and the spectra of background 
particles at time t = 2500a; In panel (b) and (e), the dash 
lines shows the corresponding quantity at t = 2470a; and 
t — 2530a;“g^ respectively (details are described in the text). 

seen some slight broadening of the energy spectra of the 
background plasma particles. Because the maximum en¬ 
ergy of the energetic particles found in the simulations 
is small, there are controversial explanations of this en¬ 
ergy spectrum broadening. Some researchers regard it as 
an acceleration [il, [s^, while others consider it as ther- 
malization [3l| . With the continuous incidents of plasma 
flows with much higher initial energy, our simulation en¬ 
ables the generation of stronger turbulent magnetic fields 
and shows irrefragable evidence for the existence of par¬ 
ticle acceleration. 


Second stage: shock acceleration 

Apparently the lepton jet propagates much faster than 
the ion-acoustic speed of the background plasma, which is 
estimated as around Va = \/ iriadEvel'^i) ~ 0.13c. Here 
7ad is the adiabatic coefficient which equals to 2 for a sim¬ 
ple particle in a two-dimensional frame. Since the system 
is highly disequilibrium, T^e is approximated with the 
average energy of background electrons, which is about 
7MeV according to the simulation. As a result, a shock 
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structure gradually forms in the background plasma as 
the interaction continues. Figure [Sfa) shows the ion den¬ 
sity distribution at t = 2b00uj~^ and its transverse av¬ 
erage is plotted in Figure (U^b) with the cases of 
earlier and later (shown by the blue and dark-red dash 
lines, respectively). One can see clearly that a shock 
wavefront propagates forward at a speed of vs — 0.253c, 
which gives a Mach number ^ = 1.95. And the density 
ratio between downstream and upstream of the shock is 
measured as around Udown/'^uv ^ 3.1 in Figure [D^b). The 
shock hydrodynamic jump conditions give 

'^down/'^u'p ~ T 1)/(Tad 1) (2) 

Vs = c[(l + 7ad7)p]/[l + 7 + ladP^] (3) 

for a well defined shock wavejs^, in which ')ad is 
the adiabatic coefficient with value 2 as discussed in the 
beginning of this paragraph, p denotes the relative mo¬ 
mentum (normalized by the rest mass multiplying c) of 
the downstream flow in the frame of the upstream (which 
in our case is the rest frame) and 7 is the related Lorentz 
factor. According to Figure O^f), we take p ^ 0.17 as 
a reasonable estimate, then 7 ^1.014. Put these values 
into the above formulas, we obtain Udownl'^up = 3.03 and 
Vs = 0.248c, which agree with the measurement in our 
simulation. 

Panels(c) and (d) of Figure [3] present the distributions 
of electromagnetic fields. One can see that at x < dOOde, 
the turbulence-generated electromagnetic fields discussed 
in the earlier part of the paper still remain strong al¬ 
though the ions are almost blown up by the lepton flow. 
The field intensity dramatically drops to some random 
fluctuations. However, at the place of the shock wave- 
front, the fields become strong again and form a peak. In¬ 
terestingly, we find that the peak positions of the electric 
and magnetic fields do not overlap with each other. The 
former locates in front of the wavefront; the latter locates 
behind the wavefront and exists in a broader range, which 
is because the sources of the fields are different. The elec¬ 
tric field is generated by the charge separation in front of 
the wave front while the source of the magnetic field is 
the instabilities. The fields move with the shock and then 
reflect the ions at the shock front, with a speed of twice 
of the shock wave speed, as one can confirm through the 
high branch in Figure [3l^f) for the longitudinal phases- 
pace. Figure [31(g) shows the Px — Py projections of the 
ion’s phasespace distribution. There is no further trans¬ 
verse acceleration comparing with t = 950cc;“g^, while the 
longitudinal acceleration is significantly enhanced by the 
shock. We also plot the energy spectra of the background 
electrons and ions in Figure |3](h). The spectrum of ions 
is not perfectly inverse-power as a whole now, but there 
are piecewise inverse-power parts at 10“^ — lO^MeV and 
10^ — lO^MeV, which forms a ’’knee” structure. Simi¬ 
lar phenomena are found in astronomical observations as 
welljs^, [ 35 I. The orange line plots the spectrum of ions 
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FIG. 4. (a) The space distribution of the background ion 

density, (b) the transverse averaged ion density and (c) x—px 
projections of the ions’ phasespace distribution at time t = 
2500a;for the initial lOOMeV lepton flow case. 


located at x > SOOde (shockfront). It shows that a quasi- 
monoenergetic peak is formed by the shock-reflected ions 
at about 200MeV, which is four times larger than the en¬ 
ergy of the injecting lepton flow. Here we would like not 
to discuss the details of shock acceleration since it is well 
studied in plenty of previous publications [ibl-llij . 

To check the robustness of the above ion accelera¬ 
tion processes, we have also performed another simu¬ 
lation with higher initial energy of the lepton flow of 
^ 100 MeV( 7 e+e- = 200 ) with other conditions un¬ 
changed. Similar phenomena are observed. In the first 
stage, the protons, which are accelerated by the Fermi 
H mechanism as well as momentum transfer from the 
lepton flow longitudinally, form a perfect inverse-power 
spectrum. The cut-off energy is slightly increased to 
^ 130MeV. However, this increase is caused only by the 
increase of longitudinal momentum, which implies that 
the Fermi H mechanism is not enhanced, i.e., the average 
speed of the ’’magnetic islands” generated by the insta¬ 
bilities is insensitive to the initial energy of the flow. On 
the other hand, a lepton flow with higher energy can pro¬ 
vide more longitudinal momentum for the downstream 
electrons, which significantly increases its temperature 
to 15.7MeV. As a result, the ion-acoustic speed increases 
to 0.18c and the measured shockwave speed increases to 
Vs = 0.293c. The ions reflected by this shock wave form 
a quasi-monoenergetic structure at about 100-400MeV. 
However, for this 7 e+,e- = 200 case, we observe that 
some ions in the downstream also have very high energy 
despite their small population, as shown in Figure |4l(c), 
which is slightly different from the 7 e+,e- = 100 case. 
This may lead to some injected acceleration into the 
shock from the pre-accelerated ions in the downstream 
region. We also present the space distribution of the ion 
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density in Figures SJ^a) andlU^b) for comparison with the 
corresponding panels in Figure [3l 


CONCLUSION 


In conclusion, we have studied the interaction of a 
high-energy lepton flow with a background normal ion- 
electron plasma via PIC simulation. A scenario of two- 
stage acceleration of the background ions is identified: In 
the first stage, the main mechanism is type-II Fermi ac¬ 
celeration found both in the transverse and longitudinal 
directions. The accelerated ions forms a globally inverse- 
power energy spectrum. As time passes, the accelera¬ 
tion gradually enters the second stage as a shock wave- 
front forms in the background plasma. The shock propa¬ 
gates in the longitudinal direction with a Mach number of 
about 2, which further accelerates the background ions to 
energy higher than that gained in the first stage. While 
the ion spectrum is no longer globally inverse-power in 
the second stage, it becomes piecewise inverse-power with 
a knee structure. And the results, which will be relevant 
with the interaction between cosmic flows with interstel¬ 
lar matters, are helpful to explain the origin of the cos¬ 
mic rays with inverse-power spectrum. The acceleration 
scenario presented in this work may be tested experimen¬ 
tally in the future as the high energy lepton flows can be 
created nowadays by high power lasers 3^38|. 
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